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Photon correlations are a cornerstone of Quantum Optics. Recent works [NJP 15 025019 V 033036 
(2013), PRA 90 052111 (2014)] have shown that by keeping track of the frequency of the photons, 
rich landscapes of correlations are revealed. Stronger correlations are usually found where the system 
emission is weak. Here, we characterize both the strength and signal of such correlations, through 
the introduction of the “frequency resolved Mandel parameter”. We study a plethora of nonlinear 
quantum systems, showing how one can substantially optimize correlations by combining parameters 
such as pumping, filtering windows and time delay. 
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I. INTRODUCTION 

The quantum theory of optical coherence developed 
by Glauber in the 60s mm revolutionized the field of 
Quantum Optics by identifying photon correlations as 
the fundamental characterization of light, instead of fre¬ 
quency 0. This is a great insight since coherence had 
been understood for centuries as a feature of monochro¬ 
maticity while it is now understood in terms of factoriz¬ 
ing correlators. This has been confirmed experimentally 
with the advance of new light sources (such as the laser 
or single-photon sources |4j) as well as progress in photo¬ 
detection. In the quantum picture, the frequency v of 
light is linked to the energy E of its constituting parti¬ 
cles through Planck constant: E = hu. The standard 
approach of photon correlations has consisted so far es¬ 
sentially in detecting photons from a light source as a 
function of time, disregarding their frequency. Experi¬ 
mentally this is achieved either with the original Hanbury 
Brown-Twiss [5] configuration or by detecting individual 
photons with a streak camera [6]. For stationary signals, 
the most important photon correlation is measured by 
the second order correlation function: 

(2) _ lim (1) 
t^oo ((afa)(t))((aJa)(t + r)) 

with a(t) the light field annihilation operator of the sys¬ 
tem under study at time t. If the corresponding spectral 
shape is singled-peak, the question of frequency correla¬ 
tions of the emitted photons may appear a moot point. 
We will see shortly that it is not. In many cases, nev¬ 
ertheless, the emission is multi-peaked and it is then 
clear that Eq. 0, which correlates photons regardless 
of which peak they originate from, is leaving some in¬ 
formation aside. It is natural to inquire what are the 
correlations of each peak in isolation, or what are the 
cross-correlations between peaks mil]. Experimentally, 
this is readily achieved by inserted filters in the arms of 
the Hanbury Brown-Twiss configurations IM3] or using 


a monochromator in a streak camera set-up [14]. The¬ 
oretically, the Glauber correlator must be upgraded to 
the so-called time and frequency resolved photon corre¬ 
lations, [H HSl H6]: 

4i,r 2 ( u, i> w 2; T ) = 

lim HI wy. UJ2F2 (t + - r)^2,r 2 (t + r)A UJlA\ M 

t^oo ((At Wljri A Wljri )(t))((At W2> r 2 A W2j r 2 )(t + r)) 

( 2 ) 

where A UJu r i {t) = f^^dtie^ UJi ~ Fi ^ 2 ^ t ~ tl ^a(ti) is the 
field detected at frequency c^, within a frequency win¬ 
dow T u at time t. We have recently developed a theory 
to compute such correlations m and introduced the con¬ 
cept of “two-photon correlation spectrum” (2PS) which, 
beyond correlating merely peaks, spans over all the pos¬ 
sible combinations of photon frequencies nans]. Land¬ 
scapes of correlations of unsuspected complexity are re¬ 
vealed as a result, which are averaged out in standard 
photon detection or remain hidden when constraining to 
particular (fixed) sets of frequencies. The 2PS enlarges 
the set of tools in multidimensional spectroscopy [2(11 [24] 
and reveals a new class of correlated emission, that can 
be useful for quantum information processing [25, 26], 
enhance squeezing m or for the study of the founda¬ 
tions of quantum mechanics [HUH]. When looking at 
the full picture, strong correlations turn out to originate 
from photons not part of the spectral peaks, since a peak 
results from a single-photon transition between two real 
states. Various such photons have weak correlations and 
even when they do, they are of a classical character. In 
contrast, collective transitions that require two photons 
to undertake the emission are strongly and non-classically 
correlated. Since they involve a virtual state whose en¬ 
ergy is not fixed, unlike for real states, they are emitted 
at other frequencies than those of the peaks mmm- 

Recently, the 2PS of a nontrivial quantum emitter 
has been experimentally observed [29], with spectacular 
agreement with the theory and positively identifying in 
a rich landscape of correlations the “leapfrog emission”, 
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i.e., between two real states separated by an intermedi¬ 
ate virtual one, as well as their violation of the Cauchy 
Schwarz inequalities. The emitter was a semiconductor 
quantum dot and the physical picture that of resonance 
fluorescence in the Mollow triplet regime m- Shortly be¬ 
fore that, a 2PS of spontaneous emission was measured 
from a polariton condensate [ 13 ], which features, how¬ 
ever, no quantum correlated emission and presents in¬ 
stead a simpler and smooth landscape alternating bunch¬ 
ing and antibunching as frequencies get similar or far 
apart, due to the fundamental Boson indistinguishabil- 
ity. These experiments confirm that the theory is sound 
and robust and that the physics of photon correlation is 
ripe to take advantage of the effects uncovered by their 
tagging with a frequency. For instance, the mere Purcell 
enhancement of leapfrog processes results in TV-photon 
emitters m- 

A central theme of frequency engineering is the inter¬ 
play between signal and correlations. Correlated emis¬ 
sion transiting by virtual states is a high-order process 
and is therefore much less frequent than direct emission. 
This brings the concern of the practical measurement of 
a 2PS, since this requires measuring coincidences from 
spectral windows where the system already emits very 
little. Mathematically, this difficulty is concealed for 
both classes of correlations, Eqs. ([l]) and Q alike, by 
the normalization (denominator) which balances the in¬ 
tensity of the coincidence emission (numerator), turning 
two vanishing numbers into a finite ratio. In this text, 
we address this problem and revisit the 2PS to take into 
account the available amount of signal. To do so, in 
Section [TT] we introduce the frequency-resolved Mandel 
parameter Qr 1 ,r 2 (^i? ^ 2 ; t) that combines both correla¬ 
tions and emission intensity. In the light of this new 
parameter, we revisit the two-photon correlation map at 
r = 0 19 for several paradigmatic examples in non¬ 
linear quantum optics built around a two-level system. 
Namely, we consider both its incoherent and coherent 
driving, the latter bringing the system into the Mollow 
triplet regime, already graced with its experimental ob¬ 
servation [29] . We also consider its coupling to a cavity to 
realize the Jaynes-Cummings (JC) physics, as well as the 
biexciton configuration found in, typically, quantum dot 
systems. These systems are briefly introduced all along 
the paper, but mainly to settle notations and we refer to 
the literature for the concepts attached to them as well 
as for their relevance to our problem. 

Even the Mandel parameter does not fully capture the 
problematic of the signal, since some correlations are so 
strong that they dominate over the scarcity of emission. 
In Section m we complement the information of the 
available signal with an estimate of the measuring time 
this supposes, defining a notion of valleys of accessible 
correlations. There is considerable freedom added by 
filtering photons when studying their correlations, and 
we explore various ways to optimize them. Various ap¬ 
proaches are illustrated for various systems, focusing on 
the JC model in Section El and the biexciton cascade 


in Section [V] In the JC case, we study the optimiza¬ 
tion with the intrinsic system parameters, namely the 
cavity-photon lifetime and the pumping rate, while in 
the biexciton case, we study the dependence on extrin¬ 
sic parameters, namely, the filters linewidth and/or time 
delay. Clearly, a comprehensive analysis could be given 
along such lines to any system of interest. The present 
text aims at illustrating such points in particular cases 
and leaves it to future works to combine them in the cases 
where they will be needed. 


II. FREQUENCY-RESOLVED MANDEL 
PARAMETER 

Mandel introduced for standard photon-correlations 
(that is, without the frequency information) a param¬ 
eter [32], now bearing his name, intended to correct for 
the previously mentioned normalization issue: the bal¬ 
ancing of two vanishing quantities that provide a finite- 
value correlation which is, for all practical purposes, not 
measurable since these quantities are those accessible to 
the experiment. The “Mandel parameter” is defined, for 
a stationary signal, as: 

Q(t) = n a {g (2) (r) - 1), (3) 

where n a = lim^oo (aJ(t)a(t)) is the steady state popu¬ 
lation of the detected mode. The offset by unity makes 
the Mandel parameter negative when the light is quan¬ 
tum (in the sense that it is sub-Poissonian and as such 
has no classical counterpart). The product by n a makes 
the normalization of coincidences to the average signal 
instead of, as previously, to uncorrelated coincidences. It 
conveys, therefore, a meaningful information on the mag¬ 
nitude of the available signal. Note that Q(0) =0 results 
either from the lack of correlated emission (^^(O) = 1) 
or from too little emission (n a —)> 0). In this way, the 
Mandel parameter really characterizes the amount of cor¬ 
related emission. 

Following the spirit of Mandel, we introduce a fre¬ 
quency resolved version: 

Qr u r 2 (^i,^2',T) = 

(S'rV K, w 2 ; r) - 1), (4) 

where = 4^,r<)(0) is now the 

steady state spectrum, which represents, physically, the 
amount of photons passing through the filter of linewidth 
Ti centered at c^. Here it must be emphasized that while 
Q(t) < 0 is a sufficient condition to establish the quan¬ 
tum character of the emission, as it corresponds to a 
Cauchy-Schwarz inequality (CSI) violation, there is not 
such a straightforward interpretation for the frequency 
resolved version that would read: 

[flf 2) ( w i> w 2,0)] 2 < g^\bJi,uji,T)g^\bj 2 ,uj2,T). (5) 
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Such violation of classical inequalities gives rise to their 
own landscape of correlations [25]. In contrast, the 
anticorrelation in frequency, which we will qualify as 
“frequency antibunching” in agreement with the liter¬ 
ature [33], only reflects anti-correlations of intensities, 
which can be, or not, linked to a quantum character of 
the emission. 

Our main theme in this text is illustrated in Fig. |TJ 
starting with the 2PS of the JC (b) under weak in¬ 
coherent pumping in the regime of spontaneous emis¬ 
sion [34} l35]. in which case its spectral lineshape is simply 
the Rabi doublet (a). The physical meaning of this cor¬ 
relation map is amply discussed in Ref. m • It is enough 
for our discussion to highlight the main phenomenology, 
namely the set of horizontal and vertical lines, that corre¬ 
spond to transitions between real states, and antidiagonal 
lines lji +CJ 2 = &2,± ~ Eo that correspond to two-photon 
“leapfrog” emission from the second manifold with levels 
at energies F? 2 ,± and the ground state at energy Eq. The 
transitions at the Rabi frequency ±Ri « Eg are all an¬ 
tibunched (blue on the figure), since they are dominated 
by the decay of one polariton from the lower manifold 
and one excitation cannot be split into two polaritons, 
while the lines at (y/2 ± 1 )g are mainly bunched (red on 
the figure), since they correspond to a cascade from the 
second manifold. The presence of such cascade correla¬ 
tions in a regime of low excitations, where the second 
manifold has a vanishing probability to be excited, illus¬ 
trates the somewhat artificial character of the 2PS. The 
problem really pertains to photon-correlations in general 
rather than to the inclusion of frequency, since they sim¬ 
ilarly predicts g( 2 )(0) = 0 regardless of the pumping in¬ 
tensity, that is to say, the system exhibits antibunch¬ 
ing of its overall emission however small is the probabil¬ 
ity for two excitations to be present simultaneously, and 
therefore for photon blockade to enforce the antibunch¬ 
ing [36] . The same holds for the harmonic oscillator at 
vanishing pumping which still generates bunched statis¬ 
tics g( 2 )(0) = 2 regardless of the probability to reach 
two excitations in the system. The paradox arises from 
the fact that in the limit where the probability of two- 
photon effects vanishes, so does the possibility to per¬ 
form a measurement, since there is no signal. Instead, 
if one considers the Mandel correlations, Eq. Q, that 
are shown in Fig. [l|c), one sees how the result makes 
more physical sense: most of the nonlinear features have 
disappeared or are considerably weakened in the regions 
where there is a strong signal (the correlations tend to die 
more slowly than the signal), and the remaining features 
are concentrated on antibunching between the peaks, as 
well as a trace of the bunching cascades. The leapfrog 
correlations are extremely strong, which is a general re¬ 
sult in all systems, while the antibunching background 
that dominates the 2PS profile has now disappeared. It 
is also worth noting how the autocorrelation of each peak 
(along the diagonal) has the butterfly shape due to indis- 
tinguishability bunching enforced by filtering [19 , while 
cross-correlation between the two peaks feature a struc- 
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FIG. 1 . (Colour online) (a) One-photon spectrum for the JC 
system at low pumping, exhibiting the Rabi doublet of strong 
coupling, (b) The corresponding 2PS obtained by span¬ 
ning gtfj (cji , UJ 2 ) over all frequencies (c) The corresponding 
frequency-resolved Mandel parameter Q la (cji, tc?2), retaining 
only the well observable features, (d) Q la (cji , 002) with a black 
(grey) mask superimposed to let appear only areas where 
N c > 10 5 coincidence are obtained in a time T exp = 12 h 
(1 h) for 7 a = 1/ns. Parameters: y a = O.lg, — O.OOlg and 
Pa 0.0l7a. 


tureless, and therefore neater, antibunching. This could 
be of interest for single-photon emitters [33]. 

Although the Mandel correlation spectrum, Fig. [ljc), 
appears more physical than the underlying 2PS, 
Fig- 0b), the latter still presents us with a more funda- 
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mental physical picture. Indeed, we have merely tamed 
down the features, not removed them, and it is useful 
to keep track of correlations even though they are out of 
reach of an actual experiment. In any case, the 2PS could 
still be measured ideally and should better be regarded as 
a theoretical limiting case. The 2PS indeed converges to 
a unique result in the limit of vanishing pumping, thereby 
defining an unambiguous correlation map, while its Man- 
del counterpart tends to zero and the relative importance 
of bunching versus antibunching areas in Fig.JIJc) depend 
on one’s choice of the pumping rate. Finally, it is worth 
mentioning that by the time of writing this text, the 2 PS 
of resonance fluorescence has already been measured in 
its entirety |29] even for a large splitting of the satel¬ 
lite peaks with spectral windows of little emission. It 
seems therefore reasonable that with the ever-increasing 
technological progress, all fundamental quantum optical 
emitters, even those with much smaller emission rates, 
will be likewise characterized. 


III. VALLEYS OF ACCESSIBLE 
CORRELATIONS 

While Qf(^i?^ 2 ) provides a physically sound picture 
of which regions of the 2 PS are the most favorable for 
observation, it also suffers from its own shortcomings. 
The arbitrary scale of Q makes it difficult to attach to it 
a quantitative figure of merit. In this Section, we further 
delineate the valleys of accessible correlations based on 
a down-to-earth estimate of the numbers of coincidences 
that can be extracted from the emission. 

Assuming no correlations, the possibility of detecting 
at least one coincidence in a time window At at frequen¬ 
cies uii, within the frequency windows T^, is given by: 


p=( 1 - e ~ ni )(l - e~ n >), ( 6 ) 

where rii = S f r i (^) 7 At represents the number of filtered 
photons from a source that emits at a rate 7 . For sim¬ 
plicity we always consider symmetric filters in this Sec¬ 
tion, T 1 = T 2 . Using that definition and assuming no 
further inefficiencies in detection, we can estimate the 
experimental time required to obtain a given number of 
random coincidences, 7V C , as follows: 

T exp = N c At/p. (7) 

With this, we plot the regions that would be resolved 
with increasing experimental time, T exp . For exam¬ 
ple, assuming y a = 1 ns -1 as a quantum dot figure of 
merit mu, we show in Fig. [ljd) the frequency-resolved 
Mandel parameter with a mask over the regions for which 
the number of coincidences is N c < 10 5 for a detec¬ 
tion time of T exp = 1 (gray) and 12 hours (black) for 
At = 1 /T. This shows how the regions with a sizable 
number of coincidences reduce to those involving at least 
one peak, as expected. Therefore, a first experimental 


confirmation of these results could be to keep one branch 
of the setup on one peak, and correlate its input with 
that of the other branch sweeping the entire spectrum. 
This should display transitions from antibunching, no 
correlation, strong bunching, no correlation again and 
a weaker antibunching in the autocorrelation, due to in- 
distinguishability bunching. For this set of parameters, 
the experiment would need to run stably for a longer 
time in order to collect the same amount of signal to 
observe also the leapfrog processes without intersecting 
with the peaks. There is a non-trivial interplay of the 
system parameters that helps/hinders the observation of 


correlations which we will be explored in Section IV 


Before moving on to the optimization, we review other 
examples of nonlinear systems explored in Ref. PI in 
the light of the frequency resolved Mandel parameter and 
the estimated time to resolve it. We start by the most 
basic system that displays a non-trivial map of correla¬ 
tions, namely, the incoherently pumped two-level system 
which we recover by setting g = 0 in the JC model. The 
one-photon spectrum of this system is a single Lorentzian 
peak with broadening T a = 'y a +P a as shown in Fig.[2|c). 
Its two-photon Mandel spectrum shows a butterfly shape 
of anticorrelation, typical of two-level systems [19] . By 
choosing T = 7 ^ and 7 a = 0.1 ns -1 , the analysis of the 
measuring time shows that a small region of frequency 
antibunching with N c > 10 6 would be observed within 
one hour, whereas most of the butterfly would be ob¬ 
served within 12 hours for the same threshold of counts. 
While it is much more binding to observe, if filtering far 
in the tail of a two-level system, one should indeed ob¬ 
serve bunching, against naive expectations. 


Next, we consider a resonant coherent driving of the 
two-level system, described by Hamiltonian Hd = Q a (cr+ 
a t). In the weak driving regime, the system has been 
recently exploited to design ultra-narrow single-photon 
sources [371 ? -[40]. In the strong-driving regime, which is 
the one that we focus on in this paper, the spectrum is 
the well known Mollow triplet El, as shown in Fig. [2](f ). 
Frequency-resolved correlations for this system have been 
theoretically investigated in the past mmm and even 
measured mmm before the concept of the 2 PS was 
put forward. But in both the theoretical and experi¬ 
mental contexts, this was at the particular frequencies 
of the three peaks. Notwithstanding, interesting corre¬ 
lations arise mainly outside the peaks m eu, at the 
cost of a weaker signal. This has been confirmed exper¬ 
imentally [29] with the full reconstruction of the Mollow 
2PS. Resonance fluorescence is indeed a system ideally 
suited to pioneer a comprehensive analysis of frequency 
photon correlations, since it is obtained in the strong¬ 
driving regime of an extremely quantum emitter, which 
implies a large emission of strongly correlated photons. 
Figure [2^e) shows that with 7 ^ = 0.1 ns -1 and T = 7 ^, 
only within 1 hour of experimental time, the regions with 
N c > 5 x 10 7 unveils all the horizontal and vertical grid 
of correlations and great part of the leapfrogs. It only 
takes 12 hours to reveal the complete two-photon Man- 
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FIG. 2 . (Colour online) Panel (a-b) [(c)]: Mandel two-photon 
spectra [one-photon spectrum] for an incoherently pumped 
two-level system with P a — O.Olyo-, — 0.1 ns -1 and T — 7 ^. 
In Panel (b) we introduce a black (grey) mask for the points 
where N c < 10 6 in a time T exp = 12 h (1 h). Panel (d-e) 
[(f)]: Mandel two-photon spectra [one-photon spectrum] for a 
coherently driven two-level system with = 570 -, 70 - = 0.1 
ns -1 and r = In Panel (e) we introduce a black (grey) 
mask for the points where iV c < 5 x 10 7 in a time T exp = 
12 h (1 h). Panel (g-h) [(i)]: Mandel two-photon spectra 
[one-photon spectrum] for an incoherently pumped biexciton 
system with P a = x — lOOyo-, 7 ^ = 0.1 ns -1 and T — 570 -. 
In Panel (h) we introduce a black (grey) mask for the points 
where iV c < 5 x 10 5 in a time T exp = 12 h (1 h). 


del spectrum. 

Finally, we consider a biexciton level scheme as de¬ 
scribed in Refs. [18, 47, 48 which is relevant in semicon¬ 
ductor quantum optics as it describes the typical level 
structure of quantum dots beyond the simplest two-level 
system picture [49H5T1 . Focusing on a single polarization, 
it consist of a three-level scheme as depicted in the inset 
of Fig.[2jh), with a ground, an excitonic (at energy ujh) 


and a biexcitonic state whose energy (cjb) differs from the 
sum of its excitonic constituents by x due to Coulomb in¬ 
teraction. The one-photon spectrum is then composed 
of two peaks (at energies ujh and ujbh = c jb ~ cj h ) 
which give rise to an interesting and rich landscape of 
two-photon correlations. The most prominent feature is 
the antidiagonal corresponding to the leapfrog between 
ground and biexciton state, uq T UO2 = ojb = — X (as¬ 
suming ujh = 0 as the reference energy), with potential 
for applications in the generation of entangled photon 
pairs by frequency filtering [18]. With 7 a = 0.1 ns -1 and 
setting the threshold at N c = 10 5 random coincidences, 
within one hour, the anticorrelation area and bunching 
of the one-photon transition peaks would be observable, 
whereas in 12 hours most of its leapfrog structure would 
be revealed as well, especially by filtering on the sides 
of the one-photon peaks. Due to both its fundamental 
importance and practical applications, we will return to 
the problem of optimizing the observation of the leapfrog 
processes by changing both the intrinsic parameters as 
well as the filtering ones in Section [V] 


IV. OPTIMIZATION OF CORRELATIONS IN 
THE JAYNES-CUMMINGS MODEL 

We now illustrate how to optimize photon correlations 
thanks to frequency filtering, in the particular case of the 
JC model. We do a qualitative analysis to avoid focusing 
the discussion on a particular set of experimental figures 
of merit. We consider the system parameters are vari¬ 
ables for the optimization and postpone to next Section 
(and to other systems) the optimization through extrinsic 
parameters, e.g., the filters and detection time. 

One parameter that can be easily modified is the in¬ 
coherent pump rate, P a . The example of the previous 
Section was chosen to be well into the linear regime, i.e., 
with a very small pumping rate, namely P a = 0.0l7 a . 
Increasing pumping has two interesting consequences for 
the observation of correlations: 

1. signal increases, 

2. the system enters the nonlinear regime. 

In Fig. [3j upper row, the effect of increasing pump¬ 
ing is shown for both the spectral shape (a-d) and the 
Mandel parameter resolved in frequency (e-h). The spec¬ 
tra let appear inner peaks between the Rabi doublet, 
corresponding to transitions from the higher manifolds. 
The corresponding Mandel 2PS also develops new fea¬ 
tures at the same time as the overall intensity of the 
correlations increases, from ^0.02 with P a = 0.0lq a to 
~ 0.5 with P a = 0.57 a (note the change of the color 
scales). In particular, the higher manifolds become vis¬ 
ible in antibunching only as they get populated, while 
they manifest themselves in bunching more clearly at 
low pumping. The Jaynes-Cummings fork provides a 
well-structured set of correlations between the peaks: the 
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FIG. 3. (Colour online) Panels (a-d): One-photon spectra for a JC system with y a — O.lg, — O.OOlg and increasing pumping 
P a as depicted in the legend. Panels (e-f): Mandel two-photon spectra for a JC system with the same parameters as in (a-d). 
Panels (i-1): One-photon spectra for a JC system with P a = O.Olg, 7 a = O.OOlg and increasing cavity decay rate, y a as depicted 
in the legend. Panels (m-p): Mandel two-photon spectra for a JC system with the same parameters as in (i-1). 


inner peaks emit bunched photons but are otherwise anti¬ 
bunched with each other, or with the remotest Rabi peak, 
and are uncorrelated with the other—closer—Rabi peak. 
When various manifolds are well populated, correlations 
are dominated by real-state transitions and virtual pro¬ 
cesses shy away in comparison. 

Another parameter that is less easily tuned but that 
determines the strong-coupling property is the cavity de¬ 
cay rate. In the bottom row of Fig.[3j the effect of increas¬ 
ing 7 a is shown, always keeping the system in the strong¬ 
coupling regime (y a < 4 g). This results in the structure 
smoothing out as well as the intensity of correlations dy¬ 
ing (note, again, the color scales). The absolute scale of 
the Mandel correlations indeed decreases by one order of 
magnitude, but partly because the cavity gets less popu¬ 
lated, and increasing pumping could compensate for that. 
For 7 a = 0.5g the leapfrog antidiagonal have disappeared 
and for 7 a = g, only the anticorrelation between the Rabi 
peaks have survived, surrounding a region of indistin- 
guishability bunching. To track more quantitatively how 
the frequency-resolved correlations evolve with the cav¬ 
ity decay rate, we show in Fig. PJa) the value of Q la for 
pairs of frequencies corresponding to filtering the Rabi 
peaks. To show that there is some difference due to in- 
distinguishability bunching in the case of autocorrelation, 
we present both the filtering for the same Rabi peak (in 
solid black) and for the two Rabi peaks (dashed red). 
For 7 a > 4g, when the system reaches the weak-coupling 
regime, frequency-filtered correlations collapse into a sin¬ 


gle curve since there are no longer polariton modes in the 

/Q\ 

system. It is instructive in this case to plot g\ a J (i?i, d=i2i) 
together with the standard g^ (dotted blue) as shown in 
the inset of Fig.[4^a). The difference in this case between 
filtering the same peak or cross-correlating them becomes 
significant. As already observed, frequency filtering the 
Rabi peaks improves antibunching as compared to non- 
filtered correlations, as we are discarding the frequency 
regions with bunched photons. This is another instance 
of how frequency filtering can be used to harness correla¬ 
tions. Note also that worsening the cavity quality factor 
betters the overall antibunching while it spoils the peaks 
antibunching. 


V. OPTIMIZATION OF CORRELATIONS IN A 
BIEXCITON CASCADE 

We now study photon correlations in a biexciton cas¬ 
cade. To link with the previous discussion on the JC, we 
show in Fig. [4^b) the dependence on the pumping rate 
for two configurations, letting here the frequency win¬ 
dow grow with the pumping as T = 2 ( 7 ^ + P a ). The two 
configurations are filtering the leapfrog transition at the 
biexciton frequency (in solid black) and the biexciton- 
exciton cascade (in dashed red). Both transitions are 
bunched, Qr > 0 , due to their two-photon cascade char¬ 
acter. However, the one that corresponds to the leapfrog 
process exhibits a clear optimal pumping intensity, at 
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FIG. 4 . (Colour online) (a) Q la (gl?i, 0J2 ) as a function of 7 a 
for a JC system with P a — O.Olya, 70 - = O.OOlg for pair of fre¬ 
quencies (Ri : Ri) (solid black), (R±,—Ri) (dashed red). In¬ 
set: g!fj (oji , UJ 2 ) for the same parameters, together with stan¬ 
dard photon correlation ^^(O) (dotted blue), (b) Qr(wi,W 2 ) 
as a function of P a for a biexciton cascade with the same 
parameters as in Fig. [ 2 ] and F = 2 ( 7 ^ + P a ), for pair a of 
frequencies as indicated in inset, (c) Qr ((^ 1 ,(^ 2 ) as a function 
of r, with the same parameter of Fig. [2] 


P a ~ whereas the one-photon transitions exhibit 

two local maxima at around P a ~ O. 27 G and P a ~ 87 ^, 
that follow the successive growth of the exciton and biex¬ 
citon populations. The low pumping regime leads to 
Qr —>• 0 due to the small emission whereas in the high 
pumping case this is because one recovers the standard 
photon correlations g^ 2 \ 0 ) ~ 1 . 


A. Asymmetric filters 

We discuss in more details how the filter linewidth af¬ 
fects the correlations. In Fig. |4^c), we show the depen¬ 
dence with the filters linewidth F for correlations of the 
leapfrog cascade (solid black) and the biexciton-exciton 



0.01 0.1 1 10 100 1000 


Tlha 

FIG. 5. (Colour online) Mandel parameter Qr u r 2 ( w bh, wh) 
when filtering the two peaks of the biexciton-exciton cascade 
as a function of the widths of the spectral windows. The 
diagonal corresponds to the typical case of identical filters. 
Stronger correlations and of a varying nature are however 
obtained for asymmetric filters. Parameters are the same as 
for Fig. [2] 


cascade (dashed red). The leapfrog correlations strongly 
depend on the filter linewidth, due to the virtual nature 
of their transitions, with an optimum value at around 
47 ^. The biexciton-exciton cascade also displays a maxi¬ 
mum r~ 7a , but its dependence is much weaker. In both 
cases, and in contrast with g^ J (coq,^) where smaller 
linewidths optimize leapfrog correlations, the compro¬ 
mise for signal requires larger filter linewidths to optimize 
the intensity of correlations. 

In the analysis so far, we have always considered sym¬ 
metric filters for the two frequencies: Ti = T 2 = T. How¬ 
ever, for a cascaded emission such as biexciton —exci¬ 
ton —>> ground state, it is worth exploring the situation 
where the filters are asymmetric, Fi 7 ^ T 2 - This is shown 
in Fig. [5j filtering the biexciton/exciton peaks for a situa¬ 
tion where they have the same broadening and intensity 
(P a = 7 a). Two areas of bunching/antibunching op¬ 
pose each other depending on the relative value of Fi vs 
T 2 , separated by a frontier of no-correlation that roughly 
correspond to the case of symmetric filters, showing the 
interest in lifting this limitation even when both spectral 
peaks are equal. Such a structure is typical of photon 
cascades. From the level structure, the natural order of 
the cascade makes it indeed likely to detect a photon first 
of frequency ujbh then of ujh- Since an ujbh [<^h\ pho¬ 
ton, filtered with Fi [T^], is the first [second] photon in 
the cascade, if T 2 > Fi the time spent by the photon in 
filter lji is larger than the one u )2 which is favouring the 
simultaneous detection of the two photons of the cascade, 
and therefore, yields a strong bunching. In the opposite 
regime, when Fi > T 2 , the ujbh photon spends less time 
in the filter preventing the simultaneous detection of the 

















0.03 


0 

- 0.01 


T~fa 


FIG. 6. (Colour online) Mandel parameter Qr(^s/2 + 
lj, cjb/2 — cj; r) when filtering the leapfrog processes (antidiag¬ 
onal) of the biexciton-exciton cascade and as a function of au¬ 
tocorrelation time. This shows the contrast between leapfrog 
(virtual) processes where the correlations are symmetric in r 
and with a fast decay, and real processes that, when in¬ 
tercepted by the filters, lead to characteristic antibunching- 
bunching transitions with a slow decay. Parameters are the 
same as in Fig. [2] except for T = 10 7 a . 


two photons of the cascade and therefore yielding strong 
antibunching in the Mandel parameter. There is an opti¬ 
mum value to observe either antibunching/bunching—as 
a rule of thumb, an order of magnitude difference—since 
ultimately the observations of correlations quenches for 
very broad/asymmetric filters. This loss of correlations is 
due, interestingly, to the overlap of the filtering windows. 


B. Delayed correlations 
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FIG. 7. (Colour online) (a) Mandel parameter 

Qr(^s/2,cJs/2; r) when filtering the leapfrog processes 
of the biexciton-exciton cascade and (b) Qt(ubh, wh\t) 
when filtering the two peaks, as a function of correlation 
time r and filters linewidth T. This highlights again the 
difference between real and virtual processes. An optimum 
value of the filters linewidth can be found for the case of real 
state transitions. Parameters are the same as for Fig. [2] 


We have also restricted our attention to r = 0, i.e., 
coincidences. However, particularly for cascaded emis¬ 
sion, it is to be expected that correlations are maximum 
at nonzero delay m • To condense the bulk of the in¬ 
formation into a single figure (Fig. [ 6 |, we consider the 
joint frequency- and time-resolved Mandel 2PS along the 
leapfrog antidiagonal of Fig. [ 2 ] (e-f), Qr(^’s/ 2 +cj, cj ^/2 — 
c j;t)) for T = lOy^. In this line lies the information of 
both the leapfrogs and the one-photon cascade. Leapfrog 
emission is maximum at uj = 0 m and is symmetric 
in r, which is typical of second-order processes where 
the photons, being virtual, have no time order. Due 
to this symmetry, the optimal delay to observe correla¬ 
tions in this case is r = 0 , and they decay with the filter 
timescale 1/T. Contrarily, the biexciton-exciton photon 
cascade, appearing at uj = {uobh, ^h}, is strongly asym¬ 
metric, as clearly shown in Fig. [ 6 ] It shows a bunch- 
ing/antibunching behaviour as expected for cascades of 
real transitions, where there is a definite temporal order 
in the emission. In this case, the optimal value to ob¬ 


serve strong bunching/antibunching is r ^ 1/T and the 
correlations ultimately decay in the intrinsic timescale 
of the system given by I/ 7 a . The different timescales 
where correlations survive between leapfrog and normal 
cascaded emissions is another consequence of their differ¬ 
ent physical origin. 


C. Combining the parameters 

Finally, after having explored separately both the de¬ 
pendence on the filter linewidth, T, and the temporal de¬ 
lay of the photons, r, one can naturally think of combin¬ 
ing them to optimize correlations cumulatively. In Fig.[7| 
we show the joint r and T dependence of the frequency- 
resolved Mandel parameter for the two more relevant 
cases: the leapfrog two-photon cascade in (a) and the 
consecutive one-photon transitions in (b). As previously 
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discussed, the temporal shape of the leapfrog cascade is 
a symmetric decay of correlations within timescale 1/r, 
as shown in the figure. As T increases, so does the decay 
rate of the correlations (the correlation time and the fil¬ 
ters linewidth are roughly in inverse proportion, as shown 
by the dotted lines which are l/(r/ 7 cr )). The r = 0 
correlation also strongly decreases and is eventually sur¬ 
rounded by antibunching oscillations in r, that corre¬ 
spond to the fast off-resonance one-photon transitions. 
For the set of parameters chosen here, the optimal corre¬ 
lations are to be found at r = 0 and for T ^ 37 a . Also as 
discussed previously, in the consecutive one-photon cas¬ 
cade in Fig. [T^b), the temporal shape exhibits a typical 
asymmetric bunching/antibunching shape. The pattern 
is fairly robust but can indeed be magnified by the ap¬ 
propriate choice of filters (we consider here symmetric 
filters for simplicity). The temporal decay occurs this 
time approximately within a time scale of 1 / 7 ^, while 
the maximum value for the correlations, both bunching 
and antibunching, is obtained at r ~ 1/T. In this case, 
the maximum is found at T ^ 37 ^ and r ^ 1/T. Note 
that correlations are optimised for filters with a width 
equal to the spectral peaks ( 37 ^ for our choice of param¬ 
eters). 

VI. SUMMARY AND CONCLUSIONS 

Summing up, filtering the photons emitted by a quan¬ 
tum source has a dramatic impact on their correlations. 


Strong correlations are often found in regions of the spec¬ 
trum where there is a weak emission, making their ex¬ 
perimental detection particularly difficult, since this im¬ 
plies coincidences of rare events. We have introduced 
a “frequency-resolved Mandel parameter” as well as a 
quantitative estimate of the time required to accumulate 
a given number of coincidences, to address this prob¬ 
lematic for several paradigmatic non-linear quantum sys¬ 
tems. We have shown the considerable flexibility opened 
by frequency-filtering, either in energy or in time, with 
possibilities to enhance correlations by varying filters 
linewidths (possibly asymmetrically), temporal windows 
of detections and system parameters (such as pumping). 
Depending on whether the correlations originate from 
real states transitions or involve virtual processes, differ¬ 
ent strategies should be adapted, corresponding to their 
intrinsically different character. With the recent experi¬ 
mental demonstration of the underlying physics [291 , the 
field of two-photon spectroscopy is now ripe to power ap¬ 
plications and optimize resources based on photon corre¬ 
lations. 
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